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Abstract. The very interesting magnetic properties of frustrated magnetic molecules are often hardly 
accessible due to the prohibitive size of the related Hilbert spaces. The finite-temperature Lanczos method 
is able to treat spin systems for Hilbert space sizes up to 10^. Here we first demonstrate for exactly solvable 
systems that the method is indeed accurate. Then we discuss the thermal properties of one of the biggest 
magnetic molecules synthesized to date, the icosidodecahedron with antiferromagnetically coupled spins of 
s = 1/2. We show how genuine quantum features such as the magnetization plateau behave as a function 
of temperature. 

PACS. 75.10.Jm Quantized spin models - 75.40.Mg Numerical simulation studies - 75.50.Xx Molecular 
magnets 



1 Introduction 

The magnetism of antiferromagnetically coupled and geo- 
metrically frustrated magnetic molecules [T1I2[31|31E1ISI[Z1[H1 
[9] is a fascinating subject due to the richness of phenomena 
that are observed pUlfTT] as well as due to the similarities 
that can be drawn towards extended spin systems such 
as the two-dimensional kagome lattice |51 [Tn[T^[13lll4l[T5] . 
But although molecules constitute finite-size spin systems, 
the investigation of their magnetic properties for instance 
in the Heisenberg model - as function of both temperature 
and magnetic field - is largely restricted if not impossible 
due to the enormous size of the underlying Hilbert spaces. 
Quantum Monte Carlo (QMC) calculations are of no help 
in this case since they suffer from the so-called negative- 
sign problem [IB l lTT l lTS] . Density Matrix Renormalization 
Group (DMRG) techniques provide another very powerful 
approximation mainly for one- dimensional spin systems 
such as chains [TOlEOj . The method delivers the relative 
ground states for orthogonal subspaces. Extensions to in- 
clude the approximate evaluation of excitations have been 
developed recently [5T]. Nevertheless, the whole method 
still works best for one-dimensional systems; applications 
to magnetic molecules are rare [22] , 

A method, which can treat medium size spin systems 
irrespective of their geometric structure, is the Lanczos 
method [53] . This method yields eigenstates with extremal 
eigenvalues in orthogonal subspaces with high accuracy 
and is thus able to deliver a magnetization curve at T = 
0. An extension towards T > is the finite-temperature 
Lanczos method (FTLM) Although it was applied 
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to several Heisenberg or Hubbard model systems, see e.g. 
[24j25, 26,27,28,29, , one must say, that this method is not 
yet very common. 

In this article we investigate whether the finite-tem- 
perature Lanczos method (FTLM) is applicable for the 
Heisenberg model describing magnetic molecules. To this 
end its accuracy is first compared for exactly solvable 
cases. Thanks to recent advances in the application of 
group theoretical methods, the energy spectra of spin sys- 
tems of unprecedented size can be evaluated numerically 
exactly [30]. Thus, antiferromagnetically coupled spin sys- 
tems with the geometric structure of the cuboctahedron 
and the icosahedron with s = 3/2 will serve as test cases; 
the Hilbert space dimension is 16,777,216 for both [301131] 

132. 

Finally, the finite temperature behavior of an antifer- 
romagnetically coupled spin system with the geometric 
structure of the icosidodecahedron, that is closely related 
to the kagome lattice, will be examined for s = 1/2 for 
the first time. Although its total Hilbert space dimension 
is 1,073,741,824, the finite-temperature Lanczos method is 
able to deliver the magnetization and the heat capacity as 
function of both temperature and applied magnetic field. 

The article is organized as follows. In Section [2] basics 
of the finite-temperature Lanczos method are repeated. 
Section [3] is devoted to the discussion of the accuracy of 
the method, and in Section[4]the method is applied to the 
icosidodecahedron. The article closes with a summary. 
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2 Reminder of the finite-temperature Lanczos 
method 

For the evaluation of thermodynamic properties in the 
canonical ensemble the exact partition function Z depend- 
ing on temperature T and magnetic field B is given by 

z(r,i3)^^(He-^^k) . (1) 

Here { 1 ) } denotes an orthonormal basis of the respec- 
tive Hilbert space. Following the ideas of Refs. the 
unknown matrix elements are approximated as 

(^|e-^^k)«^(HnM)e-^^""'VMk> , (2) 

n=l 

where | n(y) ) is the n-th Lanczos eigenvector starting 

from \u) as the initial vector of a Lanczos iteration, e^'' 
denotes the associated n-th Lanczos energy eigenvalue. 
The number of Lanczos steps is chosen as N^. In addi- 
tion, the complete and thus very large sum over all states 
1 1/ ) is replaced by a summation over a subset of R random 
vectors. Altogether this yields for the partition function 

I/— 1 n— 1 

Although this already sketches the general idea, it will 
always improve the accuracy if symmetries are taken into 
account as in the following formulation 

r ^ u=l n=l 

xe-^^Sr^'|(n(i.,r)|j.,r)r (4) 

Here F labels the irreducible representations of the cm- 
ployed symmetry group. The full Hilbert space is decom- 
posed into mutually orthogonal subspaces 'H{r). 
An observable would then be calculated as 

x{ni,y,r)\0\:^,r){iy,r\ni,y,r)) . (5) 

It was noted in Ref. [26; that this approximation of the 
observable 0{T, B) may contain large statistical fluctua- 
tions at low temperatures due to the randomness of the 
set of states {\v^r)}. It was shown that this can largely 
be cured by assuming a symmetrized version of Eq. ^. 
For our investigations this is irrelevant. 

The very positive experience is that even for large 
problems the number of random starting vectors as well as 
the number of Lanczos steps can be chosen rather small, 
e.g. R ~ 20, Nl ~ 100. The later sections will provide 
further evidence for this statement. 



It is foreseeable that the method does not work op- 
timally in very small subspaces or subspaces with large 
degeneracies of energy levels especially if the symmetry is 
not broken up into irreducible representations F. The un- 
derlying reason is given by the properties of the Lanczos 
method itself that fails to dissolve such degeneracies. The 
other case of small subspaces can be solved by including 
their energy eigenvalues and eigenstates exactly. 

Another technical issue is given by the fact that the 
chosen random vectors {i', F) should be mutually orthog- 
onal. Although one could orthonormalize the respective 
vectors, this is for practical purposes not really neces- 
sary. The reason is, that two vectors with random com- 
ponents are practically always orthogonal, because their 
scalar product is a sum over fluctuating terms that nearly 
vanishes especially in very large Hilbert spaces. 

Since Lanczos iterations consist of matrix vector mul- 
tiplications they can be parallelized by openMP directives. 
In our programs this is further accelerated by an analyt- 
ical state coding and an evaluation of matrix elements of 
the Heisenberg Hamiltonian "on the fly" [33] . 



3 Cuboctahedron and Icosahedron 




Fig. 1. Structure of the cuboctahedron (left) and the icosahe- 
dron (right). The bullets represent spin sites, the edges indicate 
interactions. 

Before employing an approximation it is necessary to 
estimate its accuracy by comparing to known exact re- 
sults. For this purpose we choose two highly frustrated 
model systems that have been treated numerically exactly 
[5II1[5T1I5^I5D] . In both systems, the cuboctahedron and the 
icosahedron (Fig.[T]), the spins are supposed to be mounted 
on the vertices of the body. All spins interact antiferro- 
magnetically with their nearest neighbors, i.e. along the 
edges of the body. The complete Hamiltonian of the spin 
system is given by the Heisenberg and the Zeeman term, 
i. e. 

^ = -2 ^ J,,£,; ■s,+gfiBBY, £f • (6) 

i<j i 

Jij is the exchange parameter between spins at sites i 
and j. The antiferromagnetic case discussed in this article 
corresponds to negative Jij. For the sake of simplicity it 
is assumed that all spins have the same spin quantum 
number si = ,S2 = • • • = stv = s = 3/2 as well as the same 
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Fig. 2. Zero-field differential susceptibility of the cuboctahe- 
dron with s = 3/2. The various curves depict the investigated 
scenarios Ri; Nl = 100. The exact dependence is given by the 
dots. 



Table 1. Employed number Ri of random starting states for 
the cuboctahedron as well as the icosahedron with s — 3/2: 
the columns provide the magnetic quantum number M, the 
dimensions of the subspaces H(M), and the Ri. "exact" means 
that this subspace is included completely and exactly. The data 
for negative A/ are given by the symmetry M -f-)- — M. 



the minimum at B — i?sat/3 for i?3 = 20. This minimum 
is related to the magnetization plateau with M — A^sat/3, 
see Refs. .34,35,36,37 . It seems that for smaller R the 
higher-lying density of states is not quite accurately re- 
produced in subspaces around M = Msat/3. 



g-factor, and that = J for nearest neighbors and zero 
otherwise. 



Since 



- 0, this (simple) symmetry is used 

for the finite-temperature Lanczos calculations. Table [T] 
shows, how the complete Hilbert space is decomposed into 
subspaces H(M) with total magnetic quantum number M. 
Besides the dimensions of those subspaces the table also 
lists four scenarios Ri, R2, -R3, and i?4, that are used for 
the realization of the FTLM. As mentioned earlier, small 
subspaces, here with M > 12, are treated exactly. 

Figure [2] displays the zero-field differential susceptibil- 
ity of the cuboctahedron with s — 3/2. One notices that 
the approximate result, that anyway deviates from the ex- 
act one only for 0.5 < kBT/\J\ < 3, quickly approaches 
the exact curve with increasing number R of initial states. 
Already for i? = 20 the approximation is practically indis- 
tinguishable from the exact one; an increase to i? = 100 
does not further improve this observable. 

Figure [3] shows the same observable, but this time as 
a function of the applied field for three low temperatures. 
Here one clearly observes that Ri — 1 leads to large de- 
viations at various fields. For R2 — 5 the deviations are 
smaller but still too big for a good approximation. The 
approximations for R3 = 20 and R4 = 100 are again very 
good for the very low temperature of fesT/j J| =0.1 which 
is due to the fact that low-lying levels which are dominant 
at this temperature are well approximated with iV^ — 100 
Lanczos steps. But for temperatures of the order of the 
exchange interaction deviations can be observed around 



The magnetic properties of the icosahedron with spin 
s = 3/2 are discussed in Ref. [30]. Analogous to the cuboc- 
tahedron its Heisenberg Hamiltonian can be diagonalized 
numerically exactly with the help of point group symme- 
tries. Here we would like to compare the exact zero-field 
heat capacity with the results of the finite-temperature 
Lanczos method, again for the scenarios listed in Table [l] 
The reason to choose the heat capacity and not the sus- 
ceptibility is given by the fact that the heat capacity 
has an unusual feature at kBT/\J\ « 0.5 which can be 
described as a small low-temperature Schottky peak. It 
stems from a bunch of low-lying degenerate and nearly 
degenerate energy levels [30 . In addition the main maxi- 
mum has a rather unusual shape compared to other mag- 
netic molecules where the main maximum is sharper and 
as a function of temperature drops off much more quickly 
towards the l/T^ behavior at high temperatures. 

As one can see in Fig.[4]an approximation with just one 
starting state (i?i = 1) for each subspace is neither able 
to reproduce the Schottky peak nor the main maximum. 
This is already much better for i?2 = 5 and practically 
almost perfect for R^ = 20. We would like to emphasize 
once more that this result is achieved with a very small 
number of states. As Table [l] shows, the low-Af subspaces 
assume a size of about 1.5 millions whereas the FTLM 
generates only R ■ Nl states in these subspaces, which for 
i?3 = 20 corresponds to just 2,000 states. The calculations 
with i?4 — 100 practically coincide with those for i?3 = 20; 
for the little Schottky peak the accuracy is even further 
improved. 
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Fig. 3. Differential susceptibility of the cuboctahedron with 
s = 3/2 for three temperatures. The various curves depict the 
investigated scenarios Ri; ~ 100. The exact dependence is 
given by the dots. 



4 Icosidodecahedron 



The icosidodecahedron of antiferromagnetically coupled 
spins is a very fascinating object, see Fig. [5] for the struc- 
ture. Chemically it is realized with spins s — 5/2 (abbr. 
Peso P]), s = 3/2 (abbr. Crgo 8 ), and s = 1/2 (abbr. 
V30 [11|38]). It belongs to the class of geometrically frus- 
trated kagome-like spin systems. Due to this close relation 
the icosidodecahedron once was termed "The kagome on 
a sphere" |35]. These molecular structures exhibit gen- 
uine properties of antiferromagnetic spin systems built of 
corner sharing triangles as there are: many singlet states 
below the first triplet state, a pronounced magnetization 
plateau with M/Mgat = 1/3, and a large magnetization 
jump to saturation [39 l l3l [34l[35ll37j . 




Fig. 4. Zero-field heat capacity of the icosahedron with s — 
3/2. The various curves depict the investigated scenarios Ri; 
Nl = 100. The exact dependence is given by the dots. 




Fig. 5. Structure of the icosidodecahedron. The bullets rep- 
resent the 30 spin sites, the edges indicate the 60 exchange 
interactions. The structure is also termed Keplerate. 
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Table 2. Employed number Ri of random starting states for 
the icosidodecahedron with s = 1/2: the columns provide the 
magnetic quantum number M, the dimensions of the subspaces 
HIM), and the Ri. "exact" means that this subspace is in- 
cluded completely and exactly. 



Although these (T = 0) properties are accessible by 
means of Lanczos diagonalization in the case of s = 1/2 
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[341135] and by means of DMRG calculations for s ~ 'i/2 
and s = 5/2 [55], the evaluation of the thermal behav- 
ior, i.e. for T > 0, seemed to be impossible due to the 
prohibitive size of the Hilbert spaces. But at least for the 
icosidodecahedron with s = 1/2 the finite-temperature 
Lanczos method could be able to deliver the temperature 
dependence of the magnetic observables. Table [2] lists the 
parameters used in our FTLM calculations. As can be de- 
duced from the large dimensions of the subspaces H(M) 
such calculations are demanding. We employed the SGI 
Altix 4700 at the German Leibniz Supercomputing Cen- 
ter using openMP parallelization with up to 510 cores as 
well as our local BULL/ScaleMP computer with 128 cores. 
To provide an estimate, a run in the subspace with M = 
and i?3 — 20 together with Nl = 100 needs about a full 
day on 510 ITANIUM II cores. 
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Fig. 6. Differential susceptibility of the icosidodecahedron 
with s — 1/2 for three temperatures. The various curves depict 
the investigated scenarios Ri; Nl ~ 100. 



Figure [6] compares the results of FTLM calculations 
with three different sets of random starting states, see Ta- 
ble [2} It is astonishing how little the observable varies 
with Ri. This means that the finite-temperature Lanc- 
zos method replaces the true spectrum very effectively by 
pseudo energy eigenvalues so that gross properties are ef- 
ficiently reproduced. 

It is important to note that the pseudo energy eigen- 
values, see top of Fig. [7] have no spectroscopic meaning in 
general. Very low-lying energy levels may nearly coincide 
with the true ones due to the rapid convergence of the 
Lanczos method for extremal eigenvalues. The vast ma- 
jority of levels - together with their weights(!) - has to 
be understood as an effective representation of the energy 
level density. To make this point clearer the bottom of 
Fig. [7] displays the low-energy part of the spectrum with 
symbols whose radii represent the weights with which they 
have to be multiplied to the Boltzmann factor in the par- 
tition function. In effect the method has some similari- 
ties with the classical Wang-Landau sampling [iDlHTlH5] . 
where one also constructs an approximate density of states 
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dron with s — 1/2 for Ri = 1 and Nl ~ 100. Bottom: Low- 
energy part of the same energy spectrum, but now the size of 
the symbols represents the weight of the corresponding state. 



consisting of discretized energy intervals and their weights 
in order to later evaluate thermal properties [43] . 

The successful determination of the temperature de- 
pendence enables us to discuss several thermal properties 
of the icosidodecahedron with s — 1/2. A key question 
is the width and thermal evolution of the magnetization 
plateau with M/Mgat = 1/3. This plateau expresses it- 
self in the differential susceptibility as a dip. Classical 
calculations for the s = 5/2 case yielded a dip that is 
much narrower than the experimental findings [5]- It was 
not evident how this feature would behave in a quantum 
calculation. Figure [s] shows both the magnetization (top 
and middle) as well as the differential susceptibility (bot- 
tom). The plateau with M/Msat = 1/3 is indicated by 
an arrow. Interestingly, the plateau as well as the dip 
disappear quickly with rising temperature. Already for 
kBT/\J\ — 0.5 they are hardly visible, and the position 
of the now much broader dip is shifted to higher fields. It 
is not yet obvious - and thus will be a matter of future 
research - how this trend transfers to quantum icosido- 
decahedra with s = 5/2 or s = 3/2 and whether it would 
be sufficient to explain the experimental findings |44| . 

Another feature, that was observed for the icosidodec- 
ahedron with s = 5/2, is the near constancy of the mag- 
netization as a function of temperature for some magnetic 
fields [IS]. Figure[9]displays the theoretical magnetizations 
for low temperatures and several magnetic field strengths. 
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but now of course for s = 1/2. There are field ranges, e.g. 
B / Bsa,t ~ 0.4, . . . , 0.8, where the magnetization varies in- 
deed very little with temperature. This is also seen in the 
middle part of Fig. [Sj where the magnetization curves for 
three temperatures virtually fall on top of each other in 
the respective field interval. At the moment it is a specula- 
tion how this behavior would change for larger spins such 
as s = 5/2. One could conjecture that the field ranges of 
the plateau as well as of the rise to saturation might shrink 
relatively and thus lead to thermally stable magnetization 
in broader field intervals. 

Finally, we would like to discuss the heat capacity of 
the icosidodecahedron with s = 1/2. Figure 10 shows the 



zero-field heat capacity that is obtained for i?3 = 20 and 
Nl = 100. Since the system possesses many singlets be- 
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Fig. 10. Zero-field heat capacity of the icosidodecahedron with 
s = 1/2 and R3 = 20; Nl = 100 (solid curve). For comparison 
the zero-field differential susceptibility is given by a dashed 
curve. 



low the first triplet one expects low-temperature features 
in the specific heat curve, that are indeed clearly visible 
(solid curve). They are absent in the zero-field differential 
susceptibility (dashed curve) that is provided for compar- 
ison. In addition the main maximum of the specific heat 
is at higher temperatures than that of the susceptibility 
which points at a higher-lying density of states that shows 
up in the heat capacity but has not much impact on the 
susceptibility. 



5 Summary and Outlook 

The finite-temperature Lanczos method enabled us to eval- 
uate the thermal properties of the antiferromagnetic spin 
icosidodecahedron with s — 1/2. The magnetic suscepti- 
bility as well as the heat capacity could be determined. A 
major result is that the magnetization plateau at M/Ms^t = 
1/3 is thermally rather unstable, i.e. it disappears above 
temperatures of /csT/j J| w 0.5. 

An important open question is how our findings change 
if the spin is increased, e.g. to s = 5/2 for the iron based 
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icosidodecahedron. Intuitively one would guess that quan- 
tum features are decreased for the more classical spin. For 
the smaller but similar cuboctahedron it could be shown 
that the number of singlets below the first triplet state 
decreases with increasing spin quantum number [3 2) . This 
would have an impact on the low-temperature features 
of heat capacity. Further investigations are necessary to 
clarify such questions which are of general nature for all 
kagome-like spin systems. 
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